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Abstract 



We investigate the influence of the geomagnetic field (GF) on the Imaging Air Cherenkov Telescope technique for two northern 
(Tenerife and San Pedro Martir) and three southern (Salta, Leoncito and Namibia (the H.E.S.S.-site)) site candidates for Cherenkov 
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Telescope Array (CTA) observatories. We use the CORSIKA and sim_telarray programs for Monte Carlo simulations of 
gamma ray showers, hadronic background and the telescope response. We focus here on gamma ray measurements in the low 
energy, sub-100 GeV, range. Therefore, we only consider the performance of arrays of several large telescopes. Neglecting the 
CsJ GF effect, we find (in agreement with previous studies) that such arrays have lower energy thresholds, and larger collection areas 

below 30 GeV, when located at higher altitudes. We point out, however, that in the considered ranges of altitudes and magnetic field 

intensities, 1800-3600 m a.s.l. and 0^0 fiT, respectively, the GF effect has a similar magnitude to this altitude effect. We provide 
the trigger-level performance parameters of the observatory affected by the GF effect, in particular the collection areas, detection 
'"i tates and the energy thresholds for all five locations, which information may be useful in the selection of sites for CTA. We also 
1^ find simple scaling of these parameters with the magnetic field strength, which can be used to assess the magnitude of the GF 
II effect for other sites; in this work we use them to estimate the performance parameters for five sites: South Africa-Beaufort West, 
O USA- Yavapai Ranch, Namibia-Lalapanzi, Chile-La Silla and India-Hanle. We roughly investigate the impact of the geophysical 
conditions on gamma/hadron separation procedures involving image shape and direction cuts. We note that the change of altitude 
^ has an opposite effect at the trigger and analysis levels, i.e. gains in triggering efficiency at higher altitudes are partially balanced by 
—I losses in the separation efficiency. In turn, a stronger GF spoils both the shape and the direction discrimination of gamma rays, thus 
its effects at the trigger and analysis levels add up resulting in a significant reduction of the observatory performance. Overall, our 
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results indicate that the local GF strength at a site can be equally important as its altitude for the low-energy performance of CTA. 
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1. Introduction 

, Imaging Air Cherenkov Telescopes (lACT) detect gamma 
rays using the Cherenkov images of their electromagnetic 
phowers developing in the atmosphere. The lACT technique 
has rapidly advanced over the last 20 years (see, e.g., a re- 
ih) and, with the current generation of lACT instm- 
bQ], it is now the most accurate and sensitive detec- 
tion technique in the very high energy gamma ray astronomy. 
The Cherenkov Telescope Array (CTA), the next generation of 
lACT detectors, is expected to improve the sensitivity of present 
observatories by an order of magnitude, coveringthe energy 
range from a few tens of GeV to hundreds of TeV |]5|1 . 

In this paper we study the Geomagnetic Field (GF) effect, 
disturbing the lACT technique, which may set an inherent 



limit on the performance of CTA, especially at low energies. 
Charged particles in atmospheric showers are deflected by the 
Earth's magnetic field, which changes the geometry of the light 
poofl (as directions of photons from e^ and e are deflected in 
opposite directions) and also leads to distortions of shower im- 
ages. This GF effect was first pointed out in |01 and later studied 
theoretically, e.g., by iHillil and lEIil for the MAGIC tele- 
scope on La Palma. On the observational ground, the effect was 
confirmed by several experiments, including the non-imaging 
Cherenkov telescope in Tien-Shan mountains [12], the Mark 
6 telescope in Narrabri il3l and recently by the Pierre Auger 
Observatory lll4ll . 

Although the GF effect is well understood, its magnitude at 
the candidate sites for CTA is rather uncertain. Previous de- 
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tailed studies were done for Mark 6 ||13|] and MAGIC UM tele- 
scopes. However the obtained results were specific to the pa- 
rameters (the telescope size, pixelization, etc.) of those instm- 



' defined as the area on the ground with nearly constant density of 
Cherenkov photons 
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ments as well as to the local GF strength and direction at their 
locations. 

The GF effect affects the lACT measurements through two 
mechanisms. First, it stretches the light pool, decreasing (on 
average) the Cherenkov photon density, thus reducing the prob- 
ability of triggering the low energy showers with the lACT ar- 
ray. In addition, the distortion of shower images spoils the 
gamma/hadron separation efficiency in the analysis of regis- 
tered data and, furthermore, may affect the reconstruction of 
the gamma ray parameters such as the estimated energy and ar- 
rival direction. In this work we thoroughly investigate the first 
effect, i.e. the influence of GF on the probability of registering 
a gamma ray. We also quantify the changes of image parame- 
ters crucial for the gamma/hadron separation and the direction 
reconstruction, which allows us to discuss trends in the quality 
of the separation and reconstruction procedures for the chang- 
ing magnitude of the GF effect. We emphasize, however, that 
our quantitative results for the triggering efficiency are more 
robust than these for the image-analysis level. The latter de- 
pend on a chosen analysis procedure, therefore, one can expect 
that the analysis-level performance will be improved when op- 
timized procedures are developed for a given site and telescope 
configuration. In contrary, the constraints on the triggering ef- 
ficiency are intrinsic and will inevitably limit the observatory 
performance. 

At the trigger level, the GF effect is most pronounced in the 
low energy regime of a telescope; for gamma ray photons with 
energies well above the threshold, the Cherenkov pools have a 
sufficiently high density of photons even in strong GF. In this 
work we rely on the CTA design concepts presented in The 
proposed configurations include 3 classes of telescopes, with 
large, medium and small size, planned to provide the optimum 
performance in different energy ranges. We focus here on the 
subarrays of several large telescopes, which are dedicated for 
high sensitivity observations below 100 GeV. Then, our results 
illustrate directly the impact of the GF effect on the energy 
threshold of the whole observatory, but qualitatively similar ef- 
fects may be expected also in other classes of telescopes around 
their threshold energies. 

Apart from the local GF at an observatory location, the lACT 
performance at low energies may be significantly affected by its 
altitude (cf. 11511 ). In this paper we consider five potentially in- 
teresting sites for CTA (cf. section 11 in Js]]). We provide the 
site-specific information, however, we aim also to derive some 
more general properties. For each site we study different direc- 
tions of observation corresponding to different strength of the 
magnetic field. Then, to disentangle effects due to the altitude 
and to the GF, for each site we analyse also the lACT perfor- 
mance with vanishing GF. The considered sites span ranges of 
altitudes and the GF strength which are sufficient to derive some 
simple scaling formulas describing the dependence of perfor- 
mance parameters on both the GF strength and altitude. Thus 
obtained scaling laws allow us to estimate the magnitude of the 
GF effect for some other sites, which were indicated after the 
completion of our computations (see Sec.|5]i. 
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Table 1: The geophysical data for Argentina-Salta (A-S), Argentina-Leoncito 
(A-L), Mexico-San Pedro Martir (M), Spain-Tenerife (S) and Namibia- 
H.E.S.S. (N) sites. The values of H, Z, S, i and Bx are obtained from the data 
given by National Geophysical Data Center at www.ngdc.noaa. gov/geomag\ 
for each site we take the GF parameters for the altitude of 10 km a.s.l. The 
positive value of the declination angle S is measured from North towards East 
directions. 



2. GF at the candidate sites 

Table [llgives the geophysical data for the sites considered in 
this paper. We use the standard parametrization of the GF (see 
lfl6ll ). B = {H,0,Z), with the x-axis pointing to the local mag- 
netic north, the y-axis pointing eastward and the z-axis oriented 
downwards. 

Charged particles in the shower observed at the zenith angle 
and the azimuthal angle cf> are deflected by the Lorentz force 
which is proportional to the component of S perpendicular to 
the observation direction: 

B± = = [h^ [cos^id) + sin^id) sin^(0)) 

-HZ sin(20) cos(0) -i- sin^(6')] " . (1) 

Note that in our definition the azimuthal angle, cf>, is measured 
with respect to the local magnetic south; Table[T]gives the mag- 
netic declination angle 6 (i.e. the angle between the local mag- 
netic north and the geographic north directions). The magnetic 
inclination angle / in Table [1] is defined as an acute angle be- 
tween B and the vertical direction. For convenience / is mea- 
sured with respect to positive and negative z-axis direction for 
northern and southern sites, respectively; we use a non-standard 
definition of / in order to make direction-dependent effects in 
the observatory performance easier to assess (thus defined / cor- 
responds to the observation angle 9). 

Note that, in general, the GF is weaker in the southern sites 
(considered here) than in the northern ones. The Argentinian 
sites are located in the region of the South Atlantic Anomaly 
(SAA). The Namibian site is located at a rather large distance 
from the minimum of SAA but still it has a rather weak GF. The 
strongest GF occurs at the Mexican site. 

Fig.[T]illustrates the dependence of the transverse GF on the 
observation direction. The change of the sign of the Z com- 
ponent between the southern and northern hemisphere results 
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Figure 1: Dependence of the transverse component of GF, B±, on the observa- 
tion direction, (a) Dependence on <f> for fixed 6=10°, 30°, 50° for the H.E.S.S. 
(the solid lines) and Tenerife (the dashed lines) sites; each line is labeled by the 
corresponding value of the Zenith angle, d. (b) The normal component of GF 
averaged over </>, (Bj_), as a function of the zenith angle for: Argentina-Salta 
(the solid line), Argentina-Leoncito (the dashed line), Namibia (the triple-dot- 
dashed line), San Pedro Martir (the dotted hne) and Tenerife (the dot-dashed 
Une). 



in an opposite dependence on the azimuthal angle in southern 
and northern sites. In the latter, the normal component of B has 
smaller values-and, hence, the GF effect is weaker-for obser- 
vation in southward directions (dashed curves in Fig. [TJ). For 
southern sites, the GF effect is weaker for observations in north- 
ward directions (solid curves in Fig.[TJ). 

Finally, we discuss the dependence on the zenith angle. As 
the absorption in the atmosphere increases with zenith an- 
gle, typically sources are not observed by lACT instruments 
with zenith angle larger than 50-60°. Large zenith angles are 
not suitable for observation at low energies, as the density of 
Cherenkov photons, for a given gamma ray energy, decreases 
significantly for 6 > 30°. For example, the energy threshold 
increases by 30% with the increase of 6 from 0° to 30°, but then 
it increases by a factor of three with the increase of 6 to 50°, 
according to simulations performed for the MAGIC I telescope 
UtIi and the HESS telescope array 1 18|. For our simulations 
we choose - 30°, which allows to study large ranges of the 
strength of the transverse GF by changing and moreover is a 
critical for observations in low energy range. 

The dependence on the zenith angle at a given site results 
from the relation between the horizontal, H, and the vertical, Z, 
intensities. At the Namibian and Mexican sites, the H compo- 
nent is larger than the Z component. At these sites, the inclina- 
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Table 2: Parameters assumed in sim.telarray for our simulations. 

tion angle of B is i 30°, see Table [T] For both Argentinian 
sites as well as for Tenerife, the horizontal intensity is much 
larger than the vertical intensity and the inclination is / ^ 60°. 

A shower developing along the magnetic field line is not 
broadened by the GF. In the Namibian and Mexican sites such 
showers are observed at = 180° and 0°, respectively, and 
X! 30°. At these sites, the change of 4> for a fixed - 30° 
yields a large amplitude of variation of B^, as shown, for the 
Namibian site, by the solid curves in Fig.[T^. At the remaining 
three sites, the amplitude of variation of B^ increases with in- 
creasing 0, with the maximum variation occurring for ^ 60°, 
as shown by the dashed curves in Fig.[T^; - 30° corresponds 
to moderate changes at these sites. 

Fig. [It shows the transverse GF averaged over (p, as a 
function of 0. As we can see, for Tenerife and the Argentinian 
sites (where H > |Z|), the average {Bj^{9)} has approximately 
the same value at all 6 < 60°. At the Namibian and Mexican 
sites (where H < |Z|), the average (5±(0)) increases signifi- 
cantly for 6 > 30°, nevertheless, (Bj^(30°)) is representative for 
zenith angles, < 30°, relevant for observations in the low en- 
ergy range. We refer to these properties in Sect. 

3. MC simulations 

We use CTA Monte Carlo tools to simulate the gamma ray, 
proton shower development and the telescope response and, 
then, we apply image-analysis procedures to the simulated data. 

The EAS (Extensive Air Shower) simulations are performed 
using the CORSIKA 6.98 code adapted for CTA flsQ; the 
code includes the description of the influence of GF on EAS. 

At each site we use the same, tropical atmospheric profile 
(see e.g. IZV]), which is appropriate for all sites considered 
in our work (the expected differences in the density of the 
Cherenkov light, resulting from differences in the air density 
profiles, do not exceed a few per cent). Concerning the atmo- 
spheric transmission, at each site we use the transmission coeffi- 
cients calculated with MODTRAN 1221] for a desert region with 
a very clear atmosphere (with the aerosol profile similar to that 
of the 'navy maritime atmosphere' discussed in li2ll ): at each 
site we use the transmission coefficients appropriate for its alti- 
tude. Obviously, the atmospheric transmission may be expected 
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to vary from one site to another, and even at a given site large 
seasonal changes may be expected, resulting in changes of the 
density of the Cherenkov light by over an order of magnitude. 
This should be an important factor for site selection, however, 
during the writing of this paper even a preliminary outcome of 
the atmospheric monitoring of the candidate sites, in terms of 
extinction, were not known (as noted in 112 ill the model used 
here reproduces very well the extinction measured on La Pakna 
during clear conditions, so it should be valid also for Tenerife). 

We simulate showers induced by primary gamma rays from 
a point-like source and by protons arriving from a cone with an 
opening full angle of 14°. The basic parameters for simulations 
of gamma rays are: 

(i) energies between 3 GeV and 1 TeV generated from a power- 
law spectrum with the photon spectral index F = 2. The 
value of r = 2 is used to achieve a sufficient number of sim- 
ulated showers in the whole energy range. However, we use 
event weights to re- weight the primary spectrum to F = 2.6 
or F = 3.5, which slopes are more typical for gamma ray 
sources. Then, most of results presented in this work corre- 
spond to F = 2.6, but we check the dependence on the value of 
F by comparing them with the results for F = 3.5. 

(ii) a fixed Zenith angle 0=30°; 

(iii) the impact points generated up to D - 800 m from the 
center of the array; we have checked that this maximum impact 
distance is sufficient to properly describe the performance of 
considered lACT arrays in the assumed range of energies. 

At each site we consider the azimuthal angles (p - 0° and 
180°, corresponding to the minimum and maximum value of 
B^. We also make a more detailed study for one specific site, 
in Namibia, by considering seven azimuthal angles between 
= 0° and 180° with uniform steps of 30°. For each site we 
consider also the case of vanishing GF, to illustrate effects due 
to the change of the altitude. The same number of 5 x 10* events 
were simulated in our analysis in each case, in particular, for 
each location and each direction of observation. 
For protons we use: 

(i) energies between 10 GeV and 3 TeV generated from a 
power-law spectrum with index F = 2 re-weighted then to 
F = 2.73; 

(ii) Zenith angle 61=30°; 

(iii) view cone of 14° (full angle); 

(iv) the impact points generated up to D = 1500 m from the 
center of the array. 

Due to small fraction of triggering events, for hadronic back- 
ground we need at least one order of magnitude more generated 
events than for gamma rays. Thus, for protons we consider 
only three sites with neglected GF, i.e. Argentina-Salta, Tener- 
ife and Namibia and one with GF (Tenerife) with angles <p = 0° 
and 180°, corresponding to the minimum and maximum value 
of B^. For each case 2x10*^ protonic events were simulated. 
For the simulations of the telescope response we use the 
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CTA sim_telarray software 11201 12311 with the default pa- 
rameter set of the telescope and camera systems assumed in the 
production-1: first CTA MC mass production (see. Fig. 18 and 
Chapt. 8 of Istl or Chapt. 6 of |24|). The production-1 param- 



Figure 2: Geometric layout of telescopes used in our simulations of gamma 
rays. The circle shows the simulated area defined by the maximum impact 
parameter. D. 



ously, the performance parameters discussed in the next section 
depend on these assumptions. However, having established a 
particular set of the telescope and camera parameters, we can 
compare the performance of an lACT observatory at various GF 
strengths or altitudes, as changes of these assumptions should 
lead to systematic improvement or worsening of the lACT per- 
formance. As an example, we have investigated various trigger 
configurations, with different number of triggered telescopes, 
A^, required to register an event. However, we found that there 
are no noticeable differences between our results for different 
A^, especially concerning the GF effect, apart from the obvious 
property of systematically higher energy thresholds and smaller 
detection rates for systems with higher A^. Then, in this paper 
we present only our results for N -2. 

We consider two telescope layouts presented in Fig. |2] 
namely arrays E(L) and K(L) with 4 and 5 Large Size Tele- 
scopes (LST), respectively. We emphasize that arrays E(L) and 
K(L), studied in this work, refer to full CTA arrays E and K 
composed of only LSTs. Full arrays E and K, which were 
commonly used in design studies for CTA, contain large num- 
ber of medium and small telescopes, additionally. Arrays E(L) 
and K(L) (see Fig.[2| are slightly offset from the center of the 
simulated area, as they are part of a full, geometrically opti- 
mised CTA layout. Nevertheless, final results are not modified 
as the simulated area is large enough to enclose all the triggered 
events. All the results presented in Sect. |4] correspond to array 
E(L); the systematic changes between arrays E(L) and K(L) are 
briefly summarized in Sect. 14.11 For the sake of simplicity we 
have assumed the same night sky backgrouncfl for each site; we 
have also assumed the same array spacing equal to 105 m (see 
Fig.|2]i, we briefly comment on this assumption in Section 4. 

We analyse the images obtained from Sim_telarray sim- 
ulations using the read_cta program (which is an internal 



eters crucial for our simulations are given in Table |2] Obvi- Hz. 



The rate of the accidental stereo events due to NSB is of the order of 100 
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Figure 3: Differential trigger rates (smoothed histograms) for Argentina-Salta 
(the solid line), Argentina-Leoncito (the dashed line), San Pedro Martir (the 
dot-dashed line), Tenerife (the dotted Hne) & Namibia (the dot-dot-dashed line), 
(a) Trigger rates without the GF; the rates differ only due to the difference of 
altitudes, (b) Trigger rates for the Northern sites with the minimum and max- 
imum GF at 6 = 30°; in upper (thicker) lines are for (p = 0° and the lower 
(thinner) lines are for (p = 1 80° . (c) Trigger rates for the Southern sites with 
the minimum and maximum GF at = 30° ; the upper (thicker) lines are for 
(p = 1 80° and the lower (thinner) lines are for ^6 = 0° . 



component of the CTA simulation package). After the im- 
age cleaning, for which we use the tail-cut^ of 5.5/11 photo- 
electrons (pe) (for cleaning algorithm details, see Sec. 4 of 
12511'). we get the image parameters discussed in Sec. 14.21 



^Pixels are considered to be a part of an image and are called boundary/core 
pixels if their signal is larger than 5.5/1 1 pe and have at least one neighbouring 
pixel with charge above 1 1 pe. 
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Figure 4: Energy threshold at the sites with neglected GF; the thresholds differ 
only due to the difference of altitudes between the sites. The upper set of points 
is for F = 2.6 and the points (from right to left) con'espond to Argentina-Salta, 
Argentina-Leoncito, San Pedro Martir, Tenerife and Namibia. The lower set of 
points is for F = 3.5 and the points are (from right to left) for Argentina-Salta, 
Tenerife and Namibia. The vertical bar in the frame represents the eiTors of 
energy threshold values (±0.5 GeV). 



4. Results 

4.1. Trigger rates, collection areas and energy thresholds 

In this section we present several trigger level parameters, 
which can be used to describe the performance of ground-based 
gamma ray detectors, cf. |26i| . 

The probability of a detection of gamma rays with energy E 
is given by the effective collection area, A(E), determined in 
MC simulations by the range of the impact parameter, D, and 
the ratio of triggered, n^, and simulated, «sin„ events, A{E) - 
[nn-{E)/nsim(E)]^D^- We calculate also the differential trigger 
rates, RiE), related toAiE) by R{E) = <i>{E)xA{E), where OCE) 
is the primary gamma ray differential energy spectrum. To set 
the normalization of R(E), we assume the primary differential 
energy spectrum of <^{E) = 2.8 x 10"'^(£'/lTeV)-- <' cm^^ s^' 



TeV ', i.e., at 1 TeV, 1% of the Crab Nebula flux (cf. ||15[1 ); 



we discuss in detail our results for such a primary spectrum 
and then we briefly compare them with the results for a softer 
spectrum, with F = 3.5. Following the conventional definition 
for lACT detectors, we determine the energy threshold, Ef^, as 
the peak energy of the differential tri gger rate0. 

To illustrate the altitude effect (cf. lll5ll27ll '). we show in Fig. 
[3^ the differential trigger rates for the sites with neglected GF. 
Below 30 GeV, the trigger rates increase with increasing alti- 
tude due to two effects. First, the atmospheric transmission for 
Cherenkov light is higher, as the air mass between the shower 
maximum and the telescope is lower, at higher altitudes. Sec- 
ond, the light pool area is smaller (an obvious geometrical ef- 
fect) and, hence, the Cherenkov photons density at distances 
close to the core axis is higher at higher altitudes. These two 
effects are important only for low energy showers, which pro- 
duce Cherenkov light intensities not exceeding significantly the 
threshold level for triggering a telescope. We can see in Fig. [3^ 
that the peak energy decreases with increasing altitude. For the 



Technically, the energy peak is determined by the mean of a Gaussian fit. 
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Figure 5: Collection area (smoothed histograms) as a function of energy for 
Argentina-Salta (the solid line), Argentina-Leoncito (the dashed line), San Pe- 
dro Martir (the dot-dashed line), Tenerife (the dotted line) & Namibia (the dot- 
dot-dashed line), (a) Collection areas at all 5 sites without the GF. (b) Collection 
areas for the Northern sites with the minimum and maximum GF at 6 = 30° ; 
in upper (thicker) lines are for = 0° and the lower (thinner) lines are for 
0=1 80° . (c) Collection areas for the Southern sites with the minimum and 
maximum GF AS = 30° ; the upper (thicker) lines are for = 1 80° and the 
lower (thinner) lines are for = 0°. 



telescope parameters and sites assumed in this work, the depen- 
dence of the energy threshold on the altitude of the array can be 
well described by a linear relation, see Fig. |4] with 



[25.4- 3.6 x(/i/lkm)] GeV, 



(2) 



where li is the altitude; superscript '0' denotes the zero mag- 
netic field. Extrapolating it to higher altitudes we get £[? - 7 



GeV at 5 km a.s.l., in approximate agreement with 111511 . The 
dependence of the collection area on altitude (see Fig. [5^) is 
significant only around the threshold and at sub-threshold en- 
ergies, e.g. at 16 GeV the collection areas at 1.8 km and 3.6 



km differ by a factor of 3, but already at 20 GeV the difference 
decreases to 30%. 

Above 30 GeV, without GF, the differential trigger rates de- 
crease, and the corresponding collection areas are smaller, with 
the increase of altitude. Primary photons with such energies 
produce a sufficiently high Cherenkov photon density to trigger 
a telescope and their detection is not much affected by atmo- 
spheric absorption or changes of photon density with altitude. 
The Cherenkov light pool is smaller at higher altitudes and, as a 
consequence, the detection probability for high energy gamma 
rays is lower. This is, however, a minor effect, e.g. at 1 TeV the 
collection areas for the largest (Salta) and smallest (Namibia) 
altitudes differ by only 20%. Moreover, this effect will not nec- 
essarily affect the full CTA array, with small size telescopes 
covering a much larger area; such a more extended array should 
be less sensitive to changes of the light-pool size with altitude. 

Overall, our results are consistent with previous studies of 
the influence of altitude on the CTA performance, see fig. 19 in 
11. 

To complete our discussion of the altitude effect, we briefly 
consider the array spacing which, after optimization for the 
height of a site, may improve the performance parameters of 
the I ACT Assuming that flie LSTs separation of 105 

m, considered in this work, is optimal for sites at the altitude 
of 2000 m a.s.l. (see Chapt. 8 in iQ] or Chapt. 6 in li24ll), and 
that the light-pool radius scales linearly with the distance to the 
shower maximum, the optimal spacing of 82 m can be estimated 
for Salta (which is the highest-altitude site considered here). We 
performed simulations for array E(L) with such a spacing at the 
Salta site and we found that the total trigger rate is larger by 
« 2%, the energy threshold is reduced by less than 1 GeV (i.e. 
the difference is comparable to the precision, ±0.5 GeV, of de- 
termining the value of Eth) and the collection area increases by 
less than 20% (and only for £ < 30 GeV) with respect to our 
results for the 105 m spacing at this site. This confirms that 
a slight improvement of the low-energy performance of lACT 
can be achieved by an optimization of the spacing and the de- 
tailed study focused on this issue should be performed after the 
site-selection decisions are made. However, the spacing effect 
is weaker than those studied in this paper and our simplified as- 
sumption of the same spacing at each site does not affect our 
overall conclusions. 

Similarly to the altitude, the GF effect reduces significantly 
the trigger rates and the collection areas, see Figs. |3j),c and|5]5,c, 
only in the low energy range; for higher energy of primary 
photons, the Cherenkov light density is sufficient even if the 
Cherenkov pool is distorted by the GF effect. For B± ^ 40 j/T, 
corresponding to (p - 180° in the northern sites, lower curves in 
Figs.[3]3 and|5}5, it has a noticeable effect up to x 100 GeV; for 
these values of B±, the collection areas at 20 GeV are reduced 
by almost a factor of 3 with respect to the corresponding val- 
ues without the GF. For the Southern sites, changes in R{E) and 
A{E) resulting from the GF effect are significant up to £ a; 60 
GeV. 

Fig. |6^ summarizes our results on the energy threshold for 
all cases calculated in our work. Fig. |6j) shows a similar sum- 
mary for the total trigger rate, Rtot, obtained by integrating /?(£■) 
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Figure 6: Energy threshold (top panel) and the total gamma trigger rates for 
a pure power-law, E^^-^ spectrum with flux at 1 TeV equal to 1 % of the Crab 
source (bottom panel), both as a function of Bx for Argentina-Salta (triangles 
& solid line), Argentina-Leoncito (crosses), San Pedro Martir (stars), Tenerife 
(circles) and Namibia (boxes). The lines with the slope from the linear fit to 
the Namibian site are shown for selected sites. The lower set of circles in panel 
(a), connected by the dot-dashed line, are for Tenerife and T = 3.5, all the 
remaining points coiTespond to F = 2.6. The vertical bar in the frame, in the 
top panel represents the errors of energy threshold values (±0.5 GeV). The 
statistical error of the total detection rates (» O.OOlphot/s) is approximately 
half a size of the markers in the bottom panel. 



over the total energy range range (from 3 GeV to 1 TeV). Both 
quantities are shown as functions of (including the cases of 
null B), rather than the azimuthal angles at the sites, as we aim 
to derive more general properties. Remarkably, both quantities 
appear to scale linearly with at a fixed altitude. By fitting a 
linear function to the data points for the Namibian site we get 
the slopes of these relations; the altitude-dependent constant 
follow from the results with neglected GF. 
For the threshold energy we get 

£th = [0.15 X {BJ\^lT) + E'^^m GeV, (3) 

where E^^ is given by Eq. (|2]l. The fitted slope gives the correct 
extrapolation up to 40 //T, see the dotted line in Fig. |6^. On 
the other hand, we note some hints for flattening of this rela- 
tion for larger h, e.g. the three data points for Salta indicate the 
slope of 0.1 GeV/yt/T (even it is lower for Leoncito); however, 
deviations from the relation determined for Namibia do not ex- 
ceed the precision, ±0.5 GeV, of determining the peak position 
in R(E). 

For the trigger rate we get 

^tot = [<o,(/0 - 0.0022 X {BJlfiT)] phot/s, (4) 



with 

Rltih) = [0.033 X (/i/lkm) -t- 0.13] phot/s. (5) 

Again, we find that this relation works well up to at least 40 yt/T. 
We also find that the same slope describes correctly all sites 
except for Salta, for which we find a slight steepening with the 
slope of 0.0027 phot/s/^t^T. 

All the above results correspond to F = 2.6. For F = 3.5 
we get a significantly lower E^h, by 4 up to 9 GeV (depending 
on h and B±), furthermore, the slopes of the linear functions 
describing the Eth(h) and Eth{B±) dependence are different than 
for F = 2.6 (see Fig.|4]and|6]l. Specifically, the dependence on 
B^ is reduced for a softer primary spectrum and for F = 3.5 
the EthiB±) function is almost flat, as shown by the dot-dashed 
line in Fig. |6^. The dependence on li is also weaker than for 
F = 2.6 and the difference between Et^ for various F decreases 
with increasing altitude (Fig. |4]i. On the other hand, the total 
trigger rates are similar for both F = 2.6 and 3.5, if we assume 
that both primary spectra have the same normalization around 
the energy threshold (i.e. at ^ 20 GeV) and, therefore, R^ot for 
F = 3.5 are not shown in Fig.|6j5 for clarity of the figure. 

The integrated trigger rates are determined much more pre- 
cisely than the peak position of the diff'erential rates. Taking 
yfn^- as a measure of the statistical error, we get the relative er- 
ror of only A/?tot/^tot ~ 4 x 10"^. Therefore, we regard the 
slopes of our linear fits determined by three data points for sites 
other than Namibia to be more reliable for the /?tot(fi±) relations 
than for £'th(Bj.). Deviations from the slopes of our fits for the 
Namibian site indicate an opposite eff'ect of large magnetic field 
at large altitudes, namely, a stronger and weaker reduction of 
the lACT performance, for /?tot and Eth, respectively. Again, 
we regard the former property to be more reliable. 

We note also that the values of R^ot and ^th for array K(L) 
(with 5 telescopes) are systematically different from those pre- 
sented in Fig. |6]for array E(L). Specifically, R^ot is larger by a 
factor of 1.2, however, the slopes of Rtot{B^) and £'th(B±) are 
the same as for array E(L). 

4.2. Image parameters 

The separation of the gamma ray signal from the dominat- 
ing hadronic background is a fundamental issue in the lACT 
technique. An eff'ective way to separate gamma rays from the 
background exploits the differences in the distributions of the 
image parameters for gamma and hadronic showers. Those dif- 
ferences are used by most of gamma/hadron separation tech- 
niques, like scaled cuts and Random Forest (e.g., Ii28i,i2^^ ). In 
this section we consider the GF and altitude effects on the two, 
most sensitive, separation parameters characterising the shape 
and orientation of shower images. 

The Hillas width and length parameters [[30l , defined as the 
second central moments calculated along the minor and ma- 
jor axes of the image, are good separation parameters, at least 
for large image size^ Namely, the mean width and length as 
functions of size (referred to as the width and length profile) are 
different for gamma rays and hadrons. 



'the size is defined as the total integrated light of the shower image 
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Figure 7: The distributions of tlie Hillas widtli (left) and lengtli (riglit) param- 
eters from our gamma ray simulations for the Tenerife site. In both figures the 
curves from bottom to top coiTespond to i? = (solid), (p = 0° (dashed) and 
(p = 180° (dotted). Multiple peaks in the width distributions appear due to 
the pixelization of the camera. The peaks' sharpness is a result of histograms 
binning. 



The broadening of images by the GF effect should be re- 
flected in the length and width distributions. Indeed, as can be 
seen in Fig. |7] the relative number of gamma ray images with 
large lengths and widths increases with increasing B±; the ef- 
fect seems to be stronger for the latter, then, we focus below on 
the width profiles. At the width values < 0.05°, their distribu- 
tion appears to be not aff'ected by the GF, which is an obvious 
(camera-dependent) effect for the camera pixel size of 0.09°. 

Fig. [8] shows the width profiles obtained from our simulated 
data by dividing the size range into bins and calculating the 
mean width, <width>, and standard deviation, cr^yidth, of the 
width distribution for each bin. For gamma rays, the increase 
of fij. at a fixed h leads to the increase of <width> at sizes 
exceeding 100 pe, e.g. <width> increases by ~ 20% with the 
increase of from to 40 fiT. The increase of <width> is 
accompanied by a similar in magnitude increase of standard de- 
viation (cr), see Fig.[8j). Neglecting the GF effect, we find that 
for gamma rays the mean width, for a given size bin, increases 
with increasing altitude, see Fig.|8ja), which effect results sim- 
ply from the decreasing distance between the telescope and the 
shower maximum, see IztIi . Note that the magnitudes of the 
changes of the gamma ray width profiles due to both effects 
(i.e. GF and altitude) are similar The changes are not dramatic, 
however, neglecting them may result in systematic inaccuracies 
e.g. in the scaled cuts technique, discussed in Section l4~4l 

Another elementary information, which must be taken into 
account for gamma/hadron discrimination, concerns the shower 
direction. In the stereo (or higher multiplicity) lACT systems, 
the arrival direction of the primary particle is estimated as the 
crossing point of the main axes of images collected from the in- 
dividual telescopes, see 13111 . In our analysis, the reconstructed 
shower direction is determined by the weighted average of all 
intersection points (i.e. for all pairs of triggered telescopes), ac- 
cording to the scheme described in ll24ll . Then, we determine 
the image parameter, defined as the angular distance between 
the simulated and the reconstructed directions of the primary 
gamma ray. 

As pointed out in previous studies, the reconstruction accu- 



0.15F^ 



I 0.1 



0.05 



0.02 



0.01 







proton 


a 


Namibia GF=0, li=1.8 km a.s.l. 

— — Tenerife GF=0, h=2.2 km a.s.l. 

Tenerife GF=MiN,h=2.2l<ma.s.i. 

Saita GF=a, h=3.6 km a.s.i. 

- - - . Tenerife GF=MAX,h=2.2 km a.s.i. 










1 


' ' " gamma-ray 

1 




Log(size/pe) 2.5 



Figure 8: The mean width for gamma rays and protons (a) and standard de- 
viation for gamma rays (b) as the function of size for Namibia without GF 
(thinner solid), Tenerife without GF (dashed), Teneiife with GF at = 0° (dot- 
dashed), Argentina-Salta without GF (thicker solid) and Tenerife with GF at 
= 180° (dotted) from bottom to top in both panels. The profiles for Argentina 
and Namibia illustrate the effect of changing the altitude only; the gamma ray 
profiles for Tenerife illustrate the effect of changing GF at a fixed altitude. 



racy is affected by both the GF strength 111 111 and the altitude 
t27.1 . The former effect involves the change of the orientation 
of the individual shower images by the GF, which obviously 
spoils the direction reconstruction. The altitude effect, stud- 
ied in detail in l27l l32ll . involves the radial distribution of the 
photon density in the light pool. Namely, the increase of the 
Cherenkov light density at the higher-altitude sites (underlying 
the enhanced trigger efficiency, studied in the previous section) 
occurs mostly close to the shower axis and, therefore, showers 
with smaller impact parameters are more frequently recorded. 
However, images of such showers have more circular shape 
and the orientation of their major axes tends to be poorly de- 
termined. 

The left panel in Fig. |9] shows the 0^ plots (the distribution 
of the squared distance between the simulated and the recon- 
structed source position) for images with Size > 100 pe; images 
with smaller sizes tend to have poorly determined orientation 
and are less suitable to illustrate effects in the reconstruction 
quality. In agreement with ill ill and 1271], our 0^ plots indicate 



that the most accurate direction reconstruction is achieved for 
the lowest altitude without GF. We can see also that the alti- 
tude effect in the 0^ plots is much weaker than the GF effect. 
E.g. the number of events reconstructed within 0^ < 0.03 deg^ 
is reduced by * 30% due to the change of Bj^ from to 40 fjT, 
and only by 10% due to the increase of altitude from 1.8 to 3.6 
km a.s.l.. This can be easily understood as the altitude affects 
the accuracy of determining the image orientation but the ac- 
tual orientation is not changed, whereas the GF does change 
the orientation. 

Obviously, the precision of the direction reconstruction de- 
termines the angular resolution, briefly discussed in Section l43] 



8 



0.02 0.04 



0.4 0.6 
[deg] 



Figure 9: 0' distribution for our MC gamma simulations (left) and the corre- 
sponding cumulative 9 (right) for Namibia without GF (thinner solid), Tenerife 
without GF (dashed), Argentina-Salta without GF (thicker solid), Tenerife with 
GF at = 0° (triple dot-dashed) and Tenerife with GF at = 180° (dotted). 
Only events with Size > 100 pe are used for both plots. 



4.3. The hadronic background 

To check how the geophysical conditions affect the hadronic- 
events detection, we apphed the analysis steps described in the 
previous sections to our MC simulations of hadronic showers 
for Salta, Namibia and Tenerife. At the trigger level, the detec- 
tion efficiency of protons is much less sensitive to changes of 
both the GF and the altitude than the detection of gamma rays. 
Without the GF, the total trigger rate of protons at Salta is only 
by 2% smaller than at the HESS-site. For Tenerife, the increase 
of Bj. from to 40 yt/T results in the decrease of the trigger 
rate of protons by 10% (as compared to the 40% decrease for 
gamma rays). 

Also the protonic images, being intrinsically broader than the 
gamma images, are much less affected by both the GF and the 
altitude. In particular, their influence on the protonic width pro- 
files appears insignificant, see Fig.|8t. 

4.4. Gamma/hadron separation 

In this section we consider the influence of the GF and al- 
titude on the two most eflicient discrimination parameters, i.e. 
the scaled width (defined below) and 0^. We investigate the 
separation quality at Salta, Namibia and Tenerife, as measured 
by the quality factor 



QF^ 



(6) 



ts' "tr 



where n^.'' and n^^'j^ are the number of gamma (hadron) events 
that passed the image cleaning and a given separation proce- 
dure, respectively. 

The scaled cuts technique uses the widths (or, similarly, 
lengths) scaled to values expected for gamma rays. In a given 
size range, the scaled width distribution is computed as Wj = 
{width - <w >)la-sf,\ith, where <w> and cr^i^jth are the mean 
width and standard deviation, as defined in Section |42] Given 
the difference in the gamma and hadron width distributions il- 
lustrated in Fig. [8] even a simple cut, e.g. accepting only events 




Figure 10: The scaled width distributions summed up over size bins, for gamma 
rays (the solid lines) and for protonic background (the dotted lines) at (a) 
Namibia without GF, (b) Salta without GF, (c) Tenerife without GF and (d) 
Tenerife with GF at = 180°. In each panel, the scaled distributions are ob- 
tained using the < h' > and tr of the gamma ray distribution for the specific 
altitude and GF. Only events with Siz,e > 150 pe are used. 



with Ws < 1, results in efficient rejection of hadrons (except for 
small sizes). 

The broadening of gamma ray images, discussed in Section 
14.21 makes them more similar to hadron images, which are neg- 
ligibly affected by both the GF and altitude. Fig. [TOl illustrates 
the related effects in Ws distributions. For gamma rays, it always 
approaches a Gaussian distribution with < w^ >— and cr - \. 
Scaling the protonic width distribution with <w> and cr„idth 
of the broadened (by GF and/or altitude) gamma rays distribu- 
tion, we get distributions of protonic Wj only weakly differing 
from that of gamma rays. This, obviously, results in a smaller 
fraction of rejected proton events. 

We quantify the related effects assuming a primary spectrum 
with F - 2.6. Fig, nn shows the values of the quality factor for 
the Ws cut, with Wj < 1 . At small sizes, the separation quality is 
poor and differences between different altitudes or GF strengths 
are small. However, the differences increase with increasing 
size (and quality factor); we note that both the increase of the 
altitude from 1.8 to 3.6 km a.s.l. and the increase of from 
to 40 yuT results in a similar reduction of the QF-as could be 
expected from magnitudes of effects illustrated in Fig.|8] 

Then, we consider a direction cut, for which we accept only 
events with (9~ < 0.05 deg^. We get a much better separation 
here, with QF > 5 even for small sizes. As pointed out in 111 ill , 
the orientation discrimination of y-rays can be significantly de- 
graded by the GF. Indeed, we note a strong reduction of the 
quality factor for the 0~ cut, e.g. from QF = 8.4 (without GF) 
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to QF ^6.9 (for = 40 //T) for size = 1000 pe at Tenerife. 
We note also that the reduction of the quality factor due to in- 
crease of from to 40 /iT is, independently of size, larger by 
a factor of a: 2.5 than the reduction due to the increase of alti- 
tude from 1 .8 to 3.6 km a.s.l., in an approximate agreement with 
the magnitude of effects illustrated in the left panel of Fig.|9] 

The increasing GF results in the reduction of both the separa- 
tion quality and the gamma ray trigger rates, therefore, one can 
expect that an overall performance should be much worse in a 
stronger magnetic field. In turn, the altitude effects act in op- 
posite directions at the trigger and analysis levels, so its effects 
may be expected to partially compensate each other The num- 
bers of both hadronic and gamma ray events simulated in this 
work are too small, by at least an order of magnitude, to perform 
a full analysis and determine the post-analysis parameters, like 
the energy-dependent sensitivity or effective area. To roughly 
measure the performance including effects at both levels, we 
define the effective quality factor similarly as in Eq. ([SJ, but us- 
ing the number of simulated gamma ray photons and hadrons. 



y,h 



n". , instead of the number of detected events, i.e. 

Sim' ' 



r 

sim 



(7) 



/«'' 



While the face value of this parameter depends on the simula- 
tion details (e.g. the maximum assumed impact parameters for 
gamma rays and protons), its relative changes track both the 
trigger and separation efficiencies; note also that it can be used 
to estimate the minimum detectable flux of gamma rays (cf. 
equation (1) in Ii33i1 ). 

We find that the effective quality factor decreases from 
QFes - 3.9 to 1 .9 with the increase of GF by 40 /^T at Tenerife, 
whereas the increase of altitude from 2.2 (Tenerife) 3.6 km a.s.l. 
(Salta) results in the increase of QF^s from 3.9 to only 4.3. This 
illustrates that the GF effect indeed strengthens with respect to 
the altitude effect in the post-analysis performance. The above 
changes of QF^s by 50% (GF) and 10% (altitude) indicate that 
the change of the altitude by 1.4km can have a similar effect 
on the overall performance as changes of the magnetic field by 
~ lOjuT, which is the typical difference between the average GF 
strengths at various sites. 

For a primary spectrum with F = 3.5 we find similar trends, 
namely QF^f^ decreases from 0.40 to 0.16 with the increase of 
GF by 40 at Tenerife whereas the increase of altitude by 
1.4 km from Tenerife to Salta results in the increase of QF^s 
from 0.40 to 0.58, i.e. the GF effect is stronger than the altitude 
effect regardless of the slope of the primary spectrum. QFeff 
is much smaller for F = 3.5 due to larger relative contribution 
of lower energy events, which are less efficiently rejected, as 
compared to F = 2.6. 

We emphasize that our analysis using QFeif is very approxi- 
mate, in particular neglecting the dependence on photon energy, 
impact parameter or size as well as neglecting information from 
other Hillas parameters or, e.g., an additional information from 
time gradient ll29ll . which could possibly improve the analysis- 
level performance. Nevertheless, we took into account the most 
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Figure 1 1 : The quality factor for tlie scaled width cut as a function of the size for 
Namibia (thin-solid), Tenerife (dashed) and Salta (thicker solid) sites without 
GF and for Tenerife at = 0° (dot-dashed) and (fi = 1 80° (dashed) with GF. 




Figure 12: Angular resolution as a function of energy for the 68% containment 
radius. Three sites are presented: Argentina-Salta without GF (thick-solid line 
with triangles), Namibia without GF (thinner solid line with boxes), Tenerife 
without GF (dashed line with empty circles), Teneiife with GF at = 0° (dot- 
dashed line with starred circles) and Tenerife with GF al (jy = 180° (dotted line 
with filled circles). The energy-dependent cuts in Size are used (e.g. 85, 300 
and 2000 pe (per telescope) for 20, 300 and 500 GeV, respectively). 



efficient discrimination parameters and the above results indi- 
cate that the GF can strongly reduce the overall CTA perfor- 
mance if appropriate correction procedures are not developed. 

4.5. Angular resolution 

The angular resolution 6>,-, is defined by the so called con- 
tainment radius, r/, containing the fraction, /, of events in the 
distribution of the reconstructed event directions, i.e. it can be 
determined by the F{0r) - f condition, where F{0) is the cu- 
mulative distribution (see right panel in Fig. |9). In Fig. [12] 
we show the angular resolution calculated for / = 0.68 with 
energy-dependent cuts in Size (e.g. 85, 300 and 2000 pe (per 
telescope) for 20, 300 and 500 GeV, respectively). The increase 
of B^ by 40 ;uT results in a very strong degradation of the angu- 
lar resolution, by a factor of three stronger than the degradation 
resulting from the increase of the altitude by 1.8 km. The in- 
fluence of the geophysical parameters on the angular resolution 
is slightly weaker at higher energies, however, the change is 
small. E.g. above the primary energy of 200 GeV the effects 
are only by a; 10% weaker than in a lowest energy regime. 
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5. Comparison of sites 

We first compare the trigger-level performance, for which we 
made the detailed simulations and which allows to compare pa- 
rameters directly related to the geophysical condition, without 
an additional dependence on the analysis procedure. The mag- 
nitude of the GF effect depends on the pointing direction of a 
telescope, therefore, for each site the performance parameters 
will change from one source to another For a rough compari- 
son of the sites, we determine the mean values, /?tot,av and ^th.av 
(given in rows ' Av' in TableO, by linear interpolation between 
the extreme values (corresponding to = 0° and - 180°) to 
the average value of the transverse GF, i.e. (Bj.(30°)). Thus de- 
rived /?tot,av and Zith.av correspond to the values of the magnetic 
field indicated by the vertical line in Fig. [TJi; as noted in Sec- 
tion 2, they give a good representation of the average GF effect 
at the considered sites. 

The southern sites are located at different altitudes with a 
similar, moderate GF and, thus, the lACT performance is pri- 
marily a function of h at these sites. Salta, located at the largest 
h, has fith.av smaller by almost ^ 7 GeV, and /?tot,av larger by 
~ 20%, than the Namibian site located at the lowest altitude. 
The two northern sites are located at similar altitudes and have 
similar {B±), as a result, differences between their Rtot.av and 
£th,av do not exceed 5%. Note also that the mean performance 
parameters in the northern sites are similar to those in Namibia. 

We emphasize, again, that parameters characterising a spe- 
cific direction of observation may significantly deviate from the 
mean values. The magnitude of these deviations is indicated by 
the extreme values of performance parameters given in rows 
'Min.' and 'Max.' in Table [3] The largest deviations occur for 
the Mexican site, due to rather large GF strength and its incli- 
nation, as noted in Section 2. 

Comparing the extreme values of 7?tot for the Mexican site 
with Rtot.av for Namibia and Salta (see Table |3]l, we find that 
the change of B± from to 40 /uT results in the decrease of 
the detection rate from 0.210 to 0.122 phot/s, and it is larger 
by ^ 50% than the decrease, from 0.198 to 0.158 phot/s, corre- 
sponding to the change of the altitude from 3.6 to 1.8 km a.s.l.. 
On the other hand, comparing changes of the threshold energy 
corresponding to such changes of B^, from 16.7 to 21.7 GeV, 
and h, from 14.4 to 21.0 GeV, we find that the change of the 
latter results in a slightly larger, by a; 2 GeV, increase of the 
threshold energy. 

Note that in the particular case of observations along the 
field lines, the Mexican site are characterised by /?tot,min = 
0.210 phot/s and Eth = 16.7 GeV similar to /?tot,av = 
0.198 phot/s and ^th.av = 14.4 GeV at Salta, indicating that 
even a moderate GF at the latter site has a noticeable effect, 
compensating the difference of altitudes of Ah -1.2 km be- 
tween these two sites. 

Using the scaling laws derived in Sec |4]we estimate in Table 
|4]the energy thresholds and trigger rates for five further sites 
which are currently considered for the location of CTA obser- 
vatories. The site in India stands out from other locations due to 
its very large altitude of 4.5 km. Neglecting the GF effect, we 
could expect Zsth ~ 9 GeV at this site. However, the GF is strong 



Parameters summary 


Site 


Northern hemisphere 


Southern hemisphere 


M 


S 


A-S 


A-L 


N 


T^tot [phot/s] (ATJtot ~ ±0.001 [phot/s]) 


Min. 


0.210 


0.172 


0.214 


0.190 


0.187 


Av. 


0.153 


0.141 


0.198 


0.170 


0.158 


Max. 


0.122 


0.123 


0.190 


0.161 


0.141 


£,h [GeV] (A£,h = ±0.5 [GeV] ) 


Min. 


16.7 


18.7 


13.7 


16.9 


19.2 


Av. 


20.0 


21.1 


14.4 


17.6 


21.0 


Max. 


21.7 


22.4 


14.7 


18.0 


22.5 



Table 3: Gamma trigger rates and the energy threshold for the lowest (Min.), 
average (specifically (Bx(30°)); Av.) and the highest (Max.) values of at 9 = 
30° for Argentina-Salta (A-S), Argentina-Leoncito (A-L), Mexico-San Pedro 
Martir (M), Spain-Tenerife (S) and Namibia-H.E.S.S. (N); A«,o, and AEn, are 
errors for i?iot and ii|h, respectively. 



at this site resulting in £'th,av ~ 14 GeV, similar to that in Salta. 
Taking into account experiments other than CTA, we note that 
a similar conclusion applies to the location of LHAASO ll34ll . 
which is planned in Tibet at a similar altitude (4.3 km) and in 
a similarly strong GF. We caution, however, that both the alti- 
tude and the GF strength at these two sites exceed the maximum 
values for which we made computations and verified our linear 
relations, so the extrapolations come with an obvious risk. 

The magnitude of the changes of the performance parame- 
ters with varying observation angle is a separate factor for eval- 
uating a site. In general, sites with H < \Z\ are characterized 
by larger changes, at least for small or moderate 9. On the 
positive side, sites with H < \Z\ have a good performance at 
small observation angles, for which is roughly equal to the 
(smaller) H component. Furthermore, the relative changes of 
performance are slightly larger at smaller altitudes (an effect 
corresponding to the decrease of the altitude-dependent factors, 
E^^{h) and ((/;), in Eqs[3]and|4| and, obviously, in a stronger 
GF. Among the considered sites, a less uniform performance 
can be expected at the northern hemisphere, where the GF is in 
general stronger Two sites with |Z| > H, in US and Mexico, are 
characterised by the largest variations, e.g. a factor of 2 changes 
of the detection rates between the southward and northward di- 
rections. The third northern site with |Z| > H, in India, also 
having the strongest GF, is located at the largest altitude which 
partially reduces the performance inhomogeneity (as explained 
above), e.g. the changes of the detection rates to a factor of 
1.5. Tenerife, with H > \Z\, has a more uniform performance 
at < 30° but stronger variations can be expected for larger 6, 
especially approaching 60°. At African sites |Z| is larger than 
H, however, the GF is weaker here than at the northern sites and 
the changes are less pronounced, e.g. variations of the detection 
rate do not exceed a factor of 1.5. The South American sites 
have |Z| < H, rather weak GF and large altitudes; as a result 
they are characterised by a very homogeneous performance. 

Concluding the trigger-level comparison, we note that 
(1) the mean parameters are similar for most of the ten con- 
sidered sites, with, in particular, isth.av around 20 GeV, and the 
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differences of /?tot,av and fith.av not exceeding 10 per cent, except 
for sites discussed in (2) below; 

(2) the Argentinian (Salta) and Indian sites provide the best per- 
spective for approaching the trigger-level energy threshold of 
10 GeV; 

(3) our scalings laws allow to estimate the amount of the per- 
formance changes for given values of H, Z and h; the largest 
performance changes characterise the US and Mexican sites, 
while the weakest changes occur at the South American sites. 

For the analysis level, we cannot discuss the energy depen- 
dent effects, in particular, the differences in the post-analysis 
energy threshold. However, the differences between the qual- 
ity factors assessed in Section 4.4 give some hints for trends 
in the post-analysis performance. In general, they indicate 
that our result at the trigger level should be taken into account 
as constraints limiting the altitude effect from above (i.e. the 
post-analysis performance should be less sensitive to altitude) 
and the GF effect from below. The differences of the average 
< Bj^ > among all considered sites do not exceed 15 fiT, and 
among sites of each hemisphere considered separately, they do 
not exceed 5 fiT. Therefore, at the trigger level the average 
performance parameters are more significantly affected by the 
altitude than the GF. When the efficiency of hadron rejection is 
taken into account, even the average GF may be similarly im- 
portant for the overall performance (very roughly estimated by 
QFeif) as the altitude. However, at the current stage any further 
statement would be premature. 

We note also that the angular resolution (assessed here with- 
out applying any separations cuts) may be markedly reduced at 
sites with stronger GF. 

6. Summary and discussion 

We have studied the influence of two geophysical factors, 
the local GF and altitude, on the low-energy performance of 
the planned CTA observatories. We quantify the changes of 
the performance parameters, which may be used to balance the 
geophysical conditions against other criteria for site selection. 
We derive an approximately linear scaling of the trigger-level 
threshold energy and detection rates with both B± and h, which 
is a novel result. 

Previous detailed studies of the GF effect were performed 
only for two sites with rather large magnetic field, namely, the 
MAGlC-site in La Palma, with the GF similar to that in Tener- 
ife, and Narrabri with even larger GF. The magnitude of the 
effect on the energy threshold in moderate GF, or its compari- 
son to the altitude effect, remained rather uncertain. Recently a 
rough estimation of the GF effect in the candidate CTA sites 
was independently done in Bachelor thesis of Maria Krause 
1 3511 . Under slightly different assumptions on configuration of 
the telescope array (with, in particular, 9 large telescopes and 80 
m spacing) she derived the difference in performance parame- 
ters between the southern and northern directions of observa- 
tion approximately consistent with ours. 

Regarding the altitude effect, the possibility of the detection 
of ~ 10 GeV photons at altitudes of ~ 5 km was considered al- 



Site 


B 


U-Y 


N-L 


C 


I 


Localization 


Lat. n 


-32.2 


35.1 


-26.6 


-29.3 


32.8 


Lon. [°] 


22.5 


-112.9 


16.6 


-70.7 


79.0 


Height 


1.6 


1.7 


1.76 


2.4 


4.5 


[km a.s.l.] 












Local magnetic field S - (H, 0, Z) [j/T] 


H 


11.0 


23.6 


11.0 


20.6 


31.7 


Z(+i) 


-24.0 


42.8 


-24.9 


-11,3 


38.7 


E 


AO = 30°) [fiT] 


^ mm 


2.5 


1.0 


2.9 


12.2 


8.1 




14.5 


28.0 


14.9 


19.7 


33.0 


g max 


21.5 


41.8 


22.0 


23.5 


46.8 


Rtot [phot/s] 


Min. 


0.18 


0.18 


0.18 


0.18 


0.26 


Av. 


0.15 


0.12 


0.16 


0.17 


0.21 


Max. 


0.13 


0.09 


0.14 


0.16 


0.18 


£th [GeV] 


Min. 


20.1 


19.5 


19.5 


18.6 


10.4 


Av. 


21.9 


23.5 


21.4 


19.8 


14.2 


Max. 


22.9 


25.6 


22.5 


20.3 


16.2 



Table 4: Trigger rates and energy thresholds obtained from Eq. and (3) 
for South Africa-Beaufort West (B), USA- Yavapai Ranch (U-Y), Namibia- 
Lalapanzi (N-L), Chile-La Silla (C) and India-Hanle (I). The values of H, Z 
and B± are obtained from the data given by National Geophysical Data Center 
at www.ngdc.noaa.gov/geomag; for each site we take the GF parameters for the 
altitude of 10 km a.s.l.. Detailed information about particular sites are available 
inChapt. 9of d. 



ready in early simulations ^ and more recently it was empha- 
sized in il5ll . Indeed, assuming that our linear scaling of Eq. (|2| 
holds at large altitudes we find that the threshold energy below 

10 GeV could be achieved by an observatory, with parameters 
assumed for CTA, located at altitudes exceeding 4.3 km in the 
absence of GF. However, the currently considered sites with 
such altitudes in India (for CTA) and in Tibet (for LHAASO) 
have rather strong GF which prevents achieving threshold ener- 
gies below 10 GeV, except for the specific case of observations 
along the field lines. 

We made a basic estimation of the distortion of image pa- 
rameters, which are crucial for the gamma/hadron separation 
and the direction reconstruction. As previously suggested by 
1 2711 . and confirmed in our calculations, the gains in trigger ef- 
ficiency with increasing altitude are significantly reduced at the 
analysis level, as the extra photons have images which are less 
suitable for the separation. We confirm also the conclusion of 

1 1 ill that the rotation of images by the GF may be the major 
obstacle in the lACT performance. We find that the maximum 
changes of the GF and altitude considered here (i.e. 40 and 
1.8 km) have similar (in magnitude) effects in the distortion of 
width profiles and in the efficiency of the scaled width cuts. On 
the other hand, the GF effect is much stronger (by a factor of 
~ 3) in the angular resolution and in the related efficiency of 
the direction cuts. We notice that while the trigger-level effects 
are important only at low energies (around the threshold), the 
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effects in the hadron rejection efficiency or in the reconstruction 
quahty are significant at all energies. 

We also made a very basic (in particular, energy- 
independent) estimation of the overall performance, which in- 
dicates that even the average GF may be equally important as 
the altitude of a site. Furthermore, a site with strong variations 
of B± for different pointing directions will be characterised by 
a very inhomogeneous performance with respect to the obser- 
vation direction. 

Our qualitative conclusions are independent of the slope of 
the primary spectrum. The threshold energy is less sensitive to 
changes of both and h for a softer spectrum. However, both 
the GF and the altitude effects on the detection rates appear 
independent of the primary spectrum slope. 

We performed high-statistics gamma ray simulations for five 
sites and we applied our scaUng laws for five additional sites. 
Among the ten considered sites, only the Argentinian and In- 
dian sites have notably better, trigger-level performance param- 
eters than the remaining sites. In the post-analysis performance 
the GF effect will be strengthened and the altitude effect weak- 
ened, however, we note again that we regard the trigger-level 
information as more basic because it does not depend on the 
assumed analysis procedure. 
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